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1 Introduction 



The scenarios that have been prevailing on the financial markets in the last decade, 
as well as the financial and insurance applications to long-lived contracts such as 
equity- linked policies (see Brennan and Schwartz [1]), suggest that the equity models 
need to take into account for stochastic and possibly positive interest rates. We 
consider therefore the problem of pricing American options when the equity asset 
follows a log-normal type diffusion process and the interest rate risk r is described 
by the Cox, Ingersoll and Ross [6] (CIR hereafter) model: 

dr{t) = k(6 - r{t))dt + a r ^fr(f)dZ r {t) } r(0) = r > 0. 
In this framework, we consider tree methods for pricing American options. 
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The existing literature for pricing American options on lattice bivariate models is 
based mainly on the technique introduced by Nelson and Ramaswamy [10], who 
show that all 1-dimensional diffusion processes can be approximated with a com- 
putationally simple binomial process by transforming the original process into a 
diffusion with unit variance. In particular, Wei [2] proposes a multidimensional 
lattice using two transformations: the first one maps both the processes (interest 
rate and equity asset) into two diffusions with unit variance and the second one 
brings to an auxiliary process which is orthogonal to the transformed interest rate 
process. This allows one to build up easily a recombining bivariate lattice to guar- 
antee the tractability of the model from a computational point of view. Hilliard, 
Schwartz and Tucker [8] propose different transformations, allowing one to treat 
with a constant variance and to eliminate the correlation between processes, and 
then proceed similarly to Wei. Nevertheless, numerical experiments show that the 
methods by Wei and by Hilliard Schwartz and Tucker are not stable and robust from 
a practical point of view whenever a T increases. 

We propose here a new lattice approach which permits accurate and efficient op- 
tion pricing numerical results without any restriction on the volatility parameter 
cr r . This procedure is structurally different from the previous ones: even if the tree 
structure is again built by using suitable transformations, the transition probabili- 
ties are independent of these maps and are set starting from the original (and not 
transformed) diffusion processes. As an immediate consequence, our approximating 
tree strongly takes into account the structure of the covariance between the spot 
rate and the asset price. Comparisons for European and American options pricing 
problems show that the proposed method is numerically very precise, efficient, sta- 
ble and robust. Moreover, the numerical efficiency of our method turns out to be 
particularly important when applied to the more elaborated numerical algorithm 
in [5], where an hybrid tree- finite difference approach is set in order to numerically 
compute American options in the Heston model. In fact, the tree component for the 
hybrid model is the one set for handling the volatility process, and turns out to be 
the tree which is presented here for the interest rate process. 

The paper is organized as follows. In Section [2J we describe the bivariate model 
for the continuous-time joint evolution of the equity and of the spot interest rate 
processes. In Section [3] we recall the algorithms introduced by Wei [H] and by 
Hilliard, Schwartz and Tucker [8]. Section @] is devoted to the description of our 
proposed tree algorithm, whose convergence is proved in Section[5j Numerical results 
and comparisons related to the evaluation tree-algorithms are the contents of Section 
El 

2 The bivariate continuous model 

We are concerned in a geometric Brownian motion describing the evolution of the 
equity value with drift driven by a square root process. So, we consider, under the 
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risk-neutral probability measure, the following dynamics for the equity value 



dS(t) 



r(t)dt + a s dZ s (t), S(0) = S > 0, 



(1) 



S(t) 



where r is the short interest rate process, as is the constant stock price volatility 
and Zs is a standard Brownian motion. The risk-neutralized process for the short 
rate is described as in the CIR model, that is 



where k is a constant representing the reversion speed, 9 is the long term reversion 
target, oy > is constant and Z r is a standard Brownian motion. When dealing with 
the CIR model, the Novikov condition should be cited: if r > and if 2k9 > of 
then a.s. the process r never hits 0. 

Finally, the Brownian noises Zs and Z r are supposed to be correlated and we let p 
denote the correlation: 



3 The Wei and Hilliard- Schwartz- Tucker proce- 
dures 

Starting from the model in Section [2} we describe here the bivariate lattices pro- 
posed by Wei [2] and by Hilliard, Schwartz and Tucker [8], following the Nelson 
and Ramaswamy [10] technique. We concentrate our attention mainly to the Wei 
procedure. Furthermore, we briefly remind the Hilliard-Schwartz- Tucker transfor- 
mations. 

In Wei p3j, both the processes S and r are firstly transformed into two diffusions 
with unit variance and then, an auxiliary process orthogonal to the transformed 
interest rate process is taken into account. This allows one to build up easily a 
recombining bivariate lattice to guarantee the tractability of the model from a com- 
putational point of view. 

The first step for the construction of a recombining bivariate tree is to transform 
processes ([T]) and into diffusions with unit variance. This is done by introducing 
the variables 



X — (log S)/as and R = 2y/r/a r 
respectively. The dynamics of X and R may be easily derived applying Ito's Lemma. 



dr(t) = k{0 - r{t))dt + a r ^/r{f)dZ r {t), r(0) = r > 0, 



(2) 



d(Zs, Z r 



)(t) = pdt. 



Hence, 



dX(t) = /ix{R{t))dt + dZ s (t), 
dR(t) = n R (R(t))dt + dZ r (t), 



X(0) = (log S ) /a s , 
R(0) = 2^/a r 



(3) 



where 



(j>x{R) 



a^/4 - al/2 



and Pr(R) 



k{A6 - R 2 al) - al 
2Ra 2 r 



(4) 
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Both the transformed processes X and R have unit variance, hence they may be 
discretized independently using recombining binomial trees. Consider now the fol- 
lowing process: 

= X(t)- P R(t) 

Standard calculations give that the processes Y(t) and R(t) have null covariance. 
In fact, the dynamics of the diffusion process Y(t) has the form 

dY{t) = fi Y {R(t))dt + dZ Y (t), Y(0) = Y = J_ [(log S ) /a s - 2p^/a r ], 

V 1 ~P 

(5) 

where 

, (m (i x {R) - PHr{R) , r x 
V>y[R) = 7= • (6) 

V 1 ~ P 

and Zy is the standard Brownian motion given by 

Zs - pZ r 



so that it is orthogonal to Z r . The underlying asset S and the spot interest rate r 
can be obtained from the state variables Y and R using the relations: 



^ = lfwL >0 S( ( )=expMVT^y W +ptf W )]. (7) 

In order to construct the discrete approximations of the processes R and Y, one 
divides the time to maturity T into N intervals of length h = T/N and, since R 
and Y are both processes with unit variance, one sets the size of each up-step and 
each down-step equal to \fh and — \fh respectively. As usual, a binomial tree is 
considered to describe the evolution of the discrete approximating processes. We 
label (0,0) the starting node where the i?-process has value R(0). After i time 
steps (i = 0, . . . , N), R may be located at one of the nodes (i, k) (k = 0, . . . , i) 
corresponding to the values 

R ik = R Q + (2k-i)Vh. (8) 

Analogously, for the discrete process approximating Y, after % time steps Y may be 
located at one of the nodes (j — 0, . . . , i) corresponding to the values 

Yij = Y + (2j-i)Vh. (9) 

Transition probabilities have to be specified to ensure the matching of the local 
drift and of the local variance between the discrete and continuous model of (Y, R). 
This will guarantee that the 2-dimensional discretized process weakly converges 
to the corresponding 2-dimensional continuous-time process given by ([3]) and (JHJ) 
(see Nelson and Ramaswamy pU] for a detailed description in the one-dimensional 
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setting). To do this one has to take into account that in some regions of the tree 
it may happen that multiple jumps are needed to satisfy properly the matching 
conditions. Hence, starting from Ri jk at time ih, the process R may jump at time 
(i + l)h to the value R i+ i^ d or Ri+i, ku , with k d and k u defined as 

k d = max{/c* : R i)k + {^ R \ k h > Ri+\,k*} and k u = k d + 1, 

in which k d = i if {k* : R i>k + (fiR)i, k h > Ri+i,k*} = and (hr)^ = /J, R (Ri,k)- One 
easily gets 

where int(x) denotes the integer part of x G R, that is, for x > then int(x) is 
the largest integer not exceeding x and for x < we set int(x) = — int(— x). Now, 
the probability that the .R-process, located at the node Ri jk , reaches Ri + i tku , is well 
defined by setting 

n (^R)i,kh + Ri : k - Ri+i,k d . 1 
Pi,k = U V A 1. 

Obviously, the probability to reach R i+l kd is 1 — p ijk - 

Concerning the F-process, one has to take into account also the behavior of R 
because here /iy = /iy(i?), see (Q. So, let us assume that the pair is located at 
(Yij,Ri^) at the time-step i. Then, for the jumps in the Y-direction, one sets 

j d = max{j* : Y id + {^ Y )i,kh > Y i+ljj *} and j u = j d + 1, 

in which j d = i if {j* : Y itj + (fi Y )i,kh > Y i+1J *} = and (fi Y )i,k = ^ Y (Ri,k)- Again, 
one has 

J, = J + mt( 

Furthermore, the transition probability that the Y-process jumps to Y i+ i j u is defined 
by setting 

- _ a V (j^rkkh + Yjj - Y i+1>jd 

Pi,j,k — U V _ v A l 

and 1 — Pi,j, k is the probability of the down-jump. 

We are now ready to describe the discrete approximation scheme for the joint evolu- 
tion of the processes Y and R by considering a bivariate tree. To this purpose, one 
takes a recombining structure for the tree by merging the two univariate binomial 
trees for the state variables R and Y. At each time step i (i = 0, . . . , N), the tree 
have (i + l) 2 nodes that we label k) corresponding to the values Ri )k and Yjj 
(k,j = 0,...,i). 

Starting from the node (i,j,k), in consideration of possible multiple jumps and 
taking into account the tree structure, the process may reach one of the following 
four nodes: 

(i + 1, j u , k u ), with probability qi, ju ,h u , 
(i + 1, j u , k d ), with probability qi, ju ,k d , 
(i + l,j d , k u ), with probability qi, jd , ku , 
(i + l,j d , k d ) , with probability qi, jd ,k d , 
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where j u , jd, k u , kd are related to multiple jumps on the tree in the Y and R direction 
respectively, and qi,j u ,k u , Qij u ,k d , %j d ,k u , Qi,j d ,k d are the associated transition probabil- 
ities. Such probabilities can be computed, due to the orthogonality of the noises 
driving the two processes, as follows 

%ju,K = Pi,kPi,j,k, qi,j u ,k d = (1 - Pi,k)Pi,j,k (10) 
%H,k u = Pi,k{ 1 ~ Pi,j,k), Qi,j d ,k d = (1 - Pi,jfc)(l - Pi,j,k)- 

The tree for the joint evolution of the processes r and S is derived simply by applying 
the maps in ([7]) to the discrete scheme just defined. Therefore, each node k) of 
the bivariate tree for (S, r) corresponds to the values 



Sij t k = exp[cr 5 (A/l - p 2 Yij + pRi,k))- (11) 

and 



otherwise. 

The resulting successor values are easily identified with rj+i^, r i+ i± d and 
Si+ij u ,k d , Si + i t j d! k u , Si + ij dt k d , and the transition probabilities are kept those defined 
in (TJl. 

Once the approximating Markov chain structure is built, the price at time of an 
American Put option with maturity T and strike K is computed by the following 
backward dynamic programming equations: 

v N,j,k = (K — SN,j,k) + 

v i}jjk = max (K - S ijjt k)+, e 



Qi,ju,k u v i+l,ju,k u "I" Qi,ju,k d Vi+l,ju,k d J t~ 

Jr( ii,j d ,ku V i+l,jd,ku + Qi,jd,k d v i+l,j d ,k d 

where v^j^, i = 0,...,n and j, k = 0,...,i, provides the American option price at 
every node (i, j, k) of the tree structure. 

Hilliard, Schwartz and Tucker [8] propose a similar procedure, which can be briefly 
recalled as follows. The first transformation (ensuring constant variance) is similar 
to the one proposed by Wei: 

X = (logS) /a s , R = 2^- 

The second transformation (removing the correlation between processes) is given by 

X 1 = a r X + R, X 2 = a r X - R. 

Then a recombining bivariate tree is constructed by means of the discrete approxima- 
tions for Xi and X 2 . The multiple jumps procedure and the transition probabilities 
are set by using a procedure which is analogous to the one just described. 
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Notice that in the procedures by Wei and by Hilliard, Schwartz and Tucker, the 
structure of the transition probabilities is of a product type, since it comes from the 
non-correlation between the processes they construct (R and Y for Wei, X\ and X2 
for Hilliard, Schwartz and Tucker). This brings to make a strongly use of the process 
R. These points, that is the product-structure for the transition probabilities and 
the use of R, make the main difference with the method we propose in this paper 
(see also next Remark [TJ . 



4 The robust tree algorithm 

In this section we introduce a new procedure for the construction of a lattice for the 
description of the approximating evolution of the pair (S, r) as in ([I])-©. Here, we 
often use the drift coefficients appearing in ([T]) and (|2J), so we recall that they are 
given by 

fis(S) — r $ an d H>r( r ) — — r) 

respectively. 

Consider first the interest rate process r. Following a remark of Tian (see [13] at 
page 100, or also |12j). the right implementation of the Nelson- Ramaswamy algo- 
rithm for the CIR process consists in setting the probabilistic structure directly on 
r, and not through R. Thus, we first construct a tree for r by transforming the 
computationally simple lattice for R, by means of the map in ffl2|) . Secondly, the 
transition probabilities are specified directly on r: the matching of the local drift 
between the discrete and the continuous model is set by using the original interest 
rate CIR process r. According to Nelson and Ramaswamy [10J, at the node (i,k), 
we set (fi r )i,k = H>r( r i,k) and we define 

k d (i, k) = max{fc* : < k* < k and r i>k + (fi r )i,kh > r i+1>k *}, (13) 
k u {i, k) = mm{k* : k+l<k*<i + l and r ijk + (fi r )i,kh < r i+1:k *} (14) 

with the understanding kd(i, k) = if {k* : < k* < k and r i)k + (fi r )i,kh > 
r i+i,k*} = and k u (i, k) = i + 1 if {k* : fc+l<A;*<i + l and r^ k + (fi r )i.kh — 
r i+ i )k *} = 0. The transition probabilities are now defined as follows: starting from 
(i, k), the probability that the process jumps to k u (i, k)) is set as 

Pi,fc = 0V A 1. (15) 

r i+l,k u (i,k) ~ r i+l,k d (i,k) 

And of course, the jump to {i + 1, kd(i, k)) happens with probability 1 — p ijk . 
Concerning the lattice for S, we proceed as follows. We first consider a computa- 
tionally simple tree-structure (Uij)ij defined as 

Uij = U + (2j-i)Vh, U = —\ogS , i = 0,...,iV, j = 0,...,i. (16) 

Then, we apply the (new) transformation 

S = e asU (17) 
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and we get a tree (Sijjy for S: Sij = e asUi ^ . By combining the two lattices, we 
obtain a 2-dimensional tree 



(Si,j, r i>k ), i = 0, . . . , N, j,k = 0,... 



i. 



It is worth to say that the transformation (IT7|) does not seem to be natural in order 
to describe the evolution of the pair. Nevertheless, this is not important. In fact, 
at this stage, we only need to set up the state-space of the Markov chain that we 
want to approximate the continuous time process (S, r), and in fact as h — > one 
really gets that the discrete state-space converges to Rl. What is important now is 
to define the transition probabilities in order to link the tree with the diffusion pair 
(5,r). 

We first set the probabilistic behavior of the jumps for S at the time-step i + 1 
given the position (Si ; j,r^ k ) by matching w.r.t. the drift fi s of S. So, we set 

(Hs)i,j,k = (J>s(Si,j,r i>k ) and we define 

jd(i,j, k) = max{j* : < j* < j and Sy + (fi s ) ijj)k h > S i+1)j *} (18) 
j u {i,j, k) = mm{j* : j + 1 < j* < i + 1 and S itj + (fi s )i,j,kh < 8 i+1 j»} (19) 

with the usual understanding jd(i,j, k) = if {j* : < j* < j and Sij + (fis)i,j,kh > 
Si+ij*} = and j u (z, j, k) = i + 1 if {j* : j + 1 < j* < i + 1 and Sy + (/is)i,j,kh < 
Si+ij«} = 0. Now, starting from (i, j, fc), the probability of an up-jump of the tree 
for S is set as 

kjk = v u-s),,iJ> + Su - SwMM A i (20) 

Si+l,j-u,(i,j,k) S i+ ij d (ij :k ) 

The down-jump obviously may happen with probability 1 — Pij jk - 
We must finally set the covariance structure for the joint evolution of the processes 
S and r on the bivariate tree. Starting from the node k), by considering possi- 
ble multiple jumps and by taking into account the tree structure, the process may 
reach one of the four nodes: (i + l,j u , k u ) with probability %j u ,k u , (i + 1, ju, kd) with 
probability qi,j u , kd , (i + l,jd, K) with probability qi )jd)ku , (i + 1, id, k d ) with probabil- 
ity %jd,k d - The transition probabilities q i)ju)k ^ qi, ju , kd , qi,j d ,k u , qi,j d ,k d are computed by 
matching (at the first order in h) the conditional mean and the conditional covari- 
ance between the continuous and the discrete processes of S and r. The matching 
conditions lead to solving the following system: 

qi,j u ,k u + qi,ju,k d = Pi,j,k 
qi,j u ,k u + qi,j d ,k u = Pi,k 

Qi,ju,k u + Qi,jd,ku + %ju,k d + Qijd,k d = 1 

m iju,k u qi,j u ,k u + m i,j u ,k d qi,j u ,k d + m ij d ,k u qi,j d ,k u + m ij d ,k d qij d ,k d = p^rsfrxk^sS^h 

(21) 

where Pij jk and pi jk are given in ( TT5|) and ( 120]) respectively and 

mi,j u ,k u = (S i+ ij u — Si t j)(r i+ i !ku — r itk ), rriij Ui k d — Wi+iju ~ Si,j)( r i+i,k d — r i,k), 

m i,j d ,k u = {Si+l,jd ~ Si,j)( r i+l,ku ~ r i,k), m i,j d ,k d = (Si+l,j d ~ Sij)(r i+1)kd — r i)k ). 

(22) 
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This completely concludes the description of the approximating process for the pair 
(S,r), and we study in Section [5] the convergence to the continuous time process. 
Now, the backward induction procedure d3J) can be applied, by using the probabilities 
solving the system (l2~Tj) . In Section E] we give numerical results from the use of this 
procedure and we compare them with the ones provided by the methods of Wei and 
of Hilliard, Schwartz and Tucker. 

Remark 1 Our procedure uses directly both the evolution of the process r and the 
correlation between S and r. Therefore, the transformed process R is exploited only 
to construct the state-space of the Markov chain approximating r, and not to build 
up the transition probabilities. This is a structural difference with what done in Wei 
IT1$ and in Hilliard, Schwartz and Tucker J2J/. Indeed, their methods are strongly 
based on the use of a bivariate diffusion process whose components are driven by 
uncorrelated noises, leading to the definition of the transition probabilities by means 
of a product. And this can be done only by handling the process R. 



5 The convergence of the robust tree 

Let us briefly recall the lattice structure for the process r: for a fixed node (i,k), 
i = 0, 1, . . . , N, k — 0, 1, . . . , i + 1, one has 

R 2 a 2 

r iik = ? ' 4 r lflt, fc >o and Ri,k = Ro + (2k-i)Vh, (23) 

with Ro = —y/r~o. We recall that kd(i, k) = max K d (i, k) and k u (i, k) = min K u (i, k), 
where 

K d (i, k) = {k* : < k* < k and r i>k + (fi r )i,kh > r i+1)k *}, (24) 
K u (i, k) = {k* : k + l<k*<i + l and r Lk + (fi r )^ k h < r i+ljk *} (25) 

with the understanding k) = if K d (i, k) = and k u (i, k) = i+1 iiK u (i, k) = 0. 
The probability that the process jumps to (i + 1, k u (i, k)) is set as 

p itk = V ^)i,kh + r itk -r i+lMiik) ^ 1 ^ 

r i+l,k u (i,k) ~ r i+l,k d (i,k) 

and of course, the jump to (i + 1, kd(i, k)) happens with probability 1 — pi )k . 
The behavior of the up and down jumps is given in the following 

Lemma 2 Let 6* < 9 and 9* > 9 be such that 

0<9,< ( 9A ^ and 0*>2(0Vr o ). 

Then there exists a positive constant h\ = hi(9*,9*,K,9,o~ r ) < 1 such that for every 
h < hi the following statements hold. 
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(i) If® — r i,k < Q*\fh then kd(i, k) = k and k u (i, k) G {k+1, . . . , Moreover, 
there exists a positive constant C* > such that 

\r i+ i,k d (i,k) ~ ri, k \ < C*h 3/4 and \r i+1Mi)k ) - r itk \ < C*h z/i . (27) 

(ii) If 6*y/h < r^ k < 9* I \fh then kd(i,k) = k and k u (i,k) = k + 1. Moreover, one 
has 

ri+i,k d (i,k)-n,k = -cr r \rr^kh+-^h and r i+1)ku ^ k) -r i;k = o r \/ r^ k h + -^h. 

(28) 

Proof. We start by considering hi G (0, 1) and in what follows, we will "calibrate" 
the value of hi. 

(i) Let us first notice that since r^ k < 9 one has (/i r )i,fc > and then for every k* < k 
one has r i+ i tk * < r^ k < r i<k + (fi r )i,kh, which gives kd(i, k) = k. Let us prove ( 127|) 
for kd{i,k). If r^jfc = then r^+i^ = as well, and (|27|) trivially holds. If instead 
< rj fc < 9*\/h, then one can have both r i+ i tk = and r i+ i :k > 0. In the first case, 
it must be 

> Ri+i )k = Ri )k — vh, 

so that one actually has < R{ ik < y/h and then r itk < a^h/A. Therefore, \r i+ i^ kd — 
r i,k\ — r i,k < of /i/4, and fl27|) holds. Consider now the second case, that is < r i<k < 
8*V~h and r i+ i tk > 0. Then it must be 

r i+ x,k = -fRi+i, k = ~f( R i,k - Vh) 2 = r iik - o r y/r iik h + 

so that 



yi+l,k,i{i,k) ^i,k\ 



a 2 . 



o r y/n^h + ^-h 



and (1271) again holds. 

Let us now discuss the up-jump. We notice that 



r i+l,i+l 



— r i,k — {^r)i,kh > tq — 9*Vh — n9h > r — 9* — n9h. 



So, by taking hi < (r — 6*)/(k6) we get K u (i, k) ^ 0, and we can proceed by looking 
for the smallest integer k* > k + 1 such that k* < i + 1 and the following statements 
hold: 

r i+ i,k* > r i>k + {fJ> r )i,kh and r i+ i ik *-i < r itk + (fi r )i,kh. (29) 

Notice that in particular the first condition gives r i+ i <k * > 0. Assume first the case 
ri+i ik *-i — 0) that is Ri+i jk * > and Ri + i )k *-i < 0. Then, 

< R i+ i, k * = Ri+i, k *-i + 2\[h < 2Vh, 
so that ri + i tk * < a^h. From the first condition in (129]) we get 

n, k < r i+ i :k * + (Hr)i,kh < (of + (fJr)i,k)K 
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so that \r i+ i^* — r itk \ < cr 2 h + (of + (fi r )i,k)h, and this proves (l27j) . So, it remains 
to study the case r i+ i tk *-i > 0. Here, we have 

u 2 u 2 
ri+i,k* = -rRi+i,k* = ~~7~ + 2V7i) 

= + 2a r y/r i+li k*-ih + a r 2 /i. 

Now, by using f l29|) we get 

< r i+1 , k * - r ijk = r i+ i ife »_i - r i)k + 2o r yjr i+ ^ k *-ih + cr^/t 



< {Hr)i,kh + 20y<i/ (r ljfc + {n r )i ik h)h + cr r /i 

and by recalling that r i fc < O^y/h, we get ( 1271) also in this last case. 

{ii) Here, we split our reasonings in two different cases: (ii.a) 8*Vh < < 6 and 
(ii.b) 9 < r i)k < 9*/Vh. 

Assume (ii.a). As for the down-jump, here we have (fJL r )i,k > 0. Since for every 
k* < k one has r i+ ^ k * < r^ k , we can immediately conclude that kd(i, k) = k. 
Concerning the up-jump, we first notice that for k* > k + 1 one has Ri + \ yk * = 
Ri,k + {2(k* — k) — \)\fh > 0, so that using the square relation with R we can write 

(2(k* — k) — I) 2 
r i+ i,fc. - r iik = (2(k* - k) - l)a r ^~ k ~h+ ^ ^ } -o 2 r h. 

Therefore, we must look for the smallest integer k* > k + 1 such that the above 
r.h.s. is larger or equal to (n r )i jk h, and this reduces to require that 

2(F - k) - 1 > 2^ + Mi > kh ~ ^ 



(J. 



So, we prove that we can choose hi such that for every h < h\ the above inequality 
holds for k* — k + 1, that is 



a 



In fact, by recalling that < (n r )i,k = — r^ k ) < k6 and r^ k > 6*y/h , we can 
write 



< ftflyg < Kg ^ 1/4 



and the statement actually holds for h < hi < (— „ 

This case can be treated similarly to the previous one, so we omit the proof. 

□ 
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Remark 3 is a consequence of the proof of Lemma [H for h < h\ we get that 
n+i,k u (i,k) — fi+i,k d (i,k) > 0, (fi r )i,kh + r i)k — rj + i jfcd (j jfe ) > and r i+ i t k u (i,k) — n,k — 
(Hr)i,kh > 0. Therefore, we can actually drop the 0V and Al appearing in &b)) and 
for h <h\ we can directly write 

(fi r )i,kh + r it k — r i+ i f k d (i,k) /„„x 
Pi,k — • l"UJ 

We also notice that whenever = one has Ti + \^ d u k \ = as well and we have 
seen that (fi r )ikh < r i+hku = r i+ i ifc „_i + 1a T ^/r i+x .k u -\h + of/i. So, since r i+hku ^i < 
n,k + (Hr)i,kh, we get 

n6h < r i+ i^ ku < (Vk9 + cr r ) 2 h 
and by substituting in l[2b}) we obtain 

K6h 



Vi.k 



r i+l,k u (i,k) 



(1 + ov/Vk 



1 



(31) 



So, we state a lower bound for the up-jump probability, which is close to 1 when the 
volatility parameter a r is is close to 0. 

The behavior of the up and down jumps for the process S are much easier to study. 
So, we recall that for a fixed node i = 0, 1, . . . , N, j — 0, 1, . . . ,2 + 1, the lattice 
structure on S is given by 

S id = e asU '<i and U itj = U + (2j - z) Vh , (32) 

with Uq = log 5*0. Here, the up and down jumps depend also on the position of 
the process r: for a fixed k — 1, . . . , i+ 1, jd(i,j, k) = max Jd(i,j, k) and j u (hj, k) = 
min J u (i,j, k), where 

Jd{i,j, k) = {j* : < j* < j and S itj + (fi s )i,j,kh > S i+1>j *}, (33) 
Ju(i,j, k) = {j* : j + 1 < j* < i + 1 and S id + (fi S )i,j,kh < S i+ ij*} (34) 

again with the understanding jd(i,j, k) — if Jd(i,j, k) — and j u (i,j, k) — i + 1 if 
J u (i,j, k) = 0. The probability that the process jumps to (z + l,j u (i,j, k)) is set as 

n w (t l s)i,j,kh + Si j - S i+l j (i j k ) 
Pi,j,k = V — A 1 (35) 

and of course, the jump to (i + l,jd(i,j, k)) happens with probability 1 — pij^- 
The behavior of the up and down jumps is given in the following 

Lemma 4 Let r* > be fixed. Then there exists /i2 = ^2(^,^5) < 1 such that for 
every h < h2 and (i, k) such that ri yk G [0, r*] one has 

j d (i,j, k) = j and j u (i,j, k) = j + 1. 

As a consequence, for h < h 2 and for every (i, k) such that ri^ k G [0, r*] one has 

S l+lM ^ k) -S id = S id (e° sVK -l) and S i+1Mi , jt k) -S i:j = S^e'^ -1). (36) 
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Proof. First of all, notice that for j* < j one has < Uij, so that Si+ij* 

Sij < < ([is)i,j,kh. This gives j d (i,j, k) = j, for every k and h, so that 



Si+i,j d (i,j,k) — S i+ ij Sij — Sij(e 



CTS 



Concerning the up-jump, it is sufficient to prove that for every h sufficiently small 
one has 

Si+lj+l — Sij > (fls)ij,kh, 

that is equivalent to 

Sij(e° sVE - 1 - r i>k h) >0. 
Since e x — 1 > x for x > 0, for r i>k < r* we can write 

Sij(e as ^ — 1 — r i>k h) > Sij(asVh — r*h) = SijVh(<Js — r*Vh), 
and the last term is positive for h < h 2 < (7 s -) 2 - Finally, 

Si+ij u (ij,k) ~~ Sij = Si-j-ij+i — Sij = Sij(e — 1). 

□ 



Remark 5 We can state a remark similar to Remark 0' in Lemma [7] we actually 
proved that for h < h 2 we have S i+1 j u (ij, k ) - S i+1 j d ^j jk ) > 0, {^ s )i,j,kh + S ijk - 
Si+ij d (ij, k ) > and S i+1 j u (ij jk ) - (/is)ij,kh - S i)k > 0. Therefore, as h < h 2 , the 
up-jump probability in can be rewritten as 

(Hs)i,j,kh + S,^ - S i+1 j d (ij tk ) 
Pi,j,k = — ~ q • i 67 ) 

°i+h3u(ij>k) °i+l,jd(i,3,k) 

We are now ready to set up and discuss the transition matrix of the bivariate Markov 
chain. We recall that, starting from the node (i,j,k), we have used the following 
notations: qij u ,k u , Qij a ,k d , Qi,j d ,k u , and qij d ,k d stand for the probability to reach (i + 
1, in, fcu), (i + 1, ju, k d ), (i + l,j d , k u ) and (z + l,j d , k d ) respectively 

Proposition 6 Let r* > and S* > be fixed and set A* = {(i,j,k) : r^ k < 
r*, Sij < S 1 *}. Let 9* be as in Lemma\E and k) £ A*. We set: 

i) if (i, j, k) £ A* and r i>k < 0*\fh then 

Qi,ju,k u ~ Pi,ju,k u Pi,k u i Qi,jv.,k d Pi,ju,k u 0- ~Pi,k u )y 

1i,jd,k u ~ 0- ~ Pi,ju,ku)Pi,k u i 1i,jd,k d ~ 0- ~ Pi,ju,k u ) ( ^ ~Pi,ku)i 

ii) if(i,j,k) £ A* and r i>k > Q^yfh then qij u>ku , qij u , kd , qij d ,k u and qij d , kd are set 
as the solutions of the linear system (12 ip . 

Then there exists hs < 1 and a positive constant C such that for every h < h% and 
(i,j,k) £ A* the above probabilities are actually well defined. 
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Proof. We fix the node (Sij,ri jk ), with (i,j,k) G A*. 

i) If r^jfc = 0, the defined transition probabilities solve system (I2ip . And if is 
small enough then the transition probabilities are supposed to be close to the ones 
in 0. So, what is said in i) is that for r^ k positive but small, we do not care the 
system and we consider the behavior in 0. 

ii) We assume here that r^ k > 6*Vh. If M denotes the 4x4 matrix underlying the 
linear system ( 12TT) . then straightforward computations give 

detM = -(Si+i,j u (i,j,k) - S i+1} j d (i t j }k ))(r i+ltku (i tk ) - r i+lthd (i,k)) 

which is non null because, by construction, both factors are positive. So, a unique 
solution i6l 4 really exists for every h. One has only to check that this actually 
gives a probability distribution, and this reduces to check that all entries of x are 
non negative. We set the solution as follows: 

= Pi,j u ,KPi,k u { 1 + 9i,jM)> %2 = Pij^kui 1 ~ Pi.fcJl 1 - QhiM), 

x 3 = (1 - Pi, ju , ku )Pi,k u (l - 9ijM)> »4 = (1 - Pi,i„,feJ(l - Pi,k u )(l + 9ijM) 

(just to be clear, we have implicity given to the four transition nodes the following 
ordering: (i + l,j u , k u ), (i + 1, j u , k d ), (i + l,j d , k u ), (i + l,j d , k d )). It is clear that x 
solves the first three equations. As for the fourth one, we take h < h± A hi, h\ and 
hi given in Lemma [2] and Lemma H] respectively, so that we can use both ( )30|) and 
( )37|) . And easy computations give 

^ _ PVr<7Sy/r\kSijh - (p r )i,k{Ps)i,j,kh 2 

l,h (A u Si,j,kPi,j,k - A d Sij :k (l - pi,j,k)) (A u r i;k p^ k - A d r i>k (l - p itk )) 

_ OreSy/?XkSijh 

(A. u Si,j,kPi,j,k - A d Sij :k (l - Pi,j, k )) (A u r itk pi :k - A d r i>k (l - p itk )) 

(Pr)i,k(Ps)i,j,kh 2 

Vr°sJrXkSijh 



X 



x \P 



in which we have set A u S i:jik = S i+1 j u ^j :k ) — S it j, A d S iij>k = S i+1: j d ^j :k ) — S i:j , 
A u r^ k = r i+1)fcu ( i)fc ) - r i)k and A d r ijk = r i+1)kd{i)k) - r iik . So, we only need to show 
that for small values of h one gets 

sup \g(p)\ < 1. 

(i,j,k)eA„ and r iik >6*\/h. 

We write 

/ \ 1 ( (Pr)i.k(Ps)i.j,kh 2 \ 

9i,j,k{p) = [p g , ■ 

Oii,j,k V Vr&SyPXkbi,jh / 

with 

(A u S i}jjk pi t j tk - A d S i jjtk (l - Pi,j,k)) {A u r itk p itk - A d r i}k (l - p i>k )) 



CTrCTsWrilSijh 
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We use (1281) . (|36|) and we write, for h small, 

AA-fc = ^e CT ^(l - e"^) > 5^(1 - e~° sVK ) = -A d Si 

a 2 ^2 

A n r ijfc = a r \fr~Ji + -j- /i > a r \fr itk h - -j- h = -A d r^ k 



so that 



&r&S\/rilSi,jh 



By recalling that 1 — e crs ^ 1 > asVh — ^f- /ie°" s , we get 



a i ,j,k > I 1 



^vWi-^ -Wi 



4 \/ r, 



in which we have used that > O^y/h. Setting c* 

tti^fc >(l - c* /i^) 2 







<7 r 


2 




~ 4^ 


/x 2 


, we get 



Therefore, we obtain 



\9i,jM\ < 



1x2 



< 



< 



(l — c* 
1 



P ~ 



(/J'r)i,k(lJ's)i,j,kh^ 



{Pr)i,k(Ps)i,j,kh 2 



(l - C* hi) 
1 



(l - c* hi) 



— 2 (\p\ + K[ / ^ ) h 1/4 + -h^ 



and since \p\ < 1, the last r.h.s. can be set less than 1 for h small enough. 
□ 

We can now state the main result. In order to do this, we set (5^ , r\ )i=o,i,...,jv the 
Markov chain running on the bivariate lattice structure, that is: 

• Sq = S and r# = r ; 

• at time ih, the state-space for the pair (S^,r^) is given by {(Sij,r i)k ) : j,k = 
0,1,. ..,i}; 

• from time i/i to time (z + l)/i the transition law on M 2 is given by 

n h (r i)fc , da) =g i j li ,fc ll 5{(5 l+1 , JU ,r l+1 , fe j}(da;) + ?ij u ,fe/{(5 i+1 , 3U) r- i+lifed )}(o?x) + 
+ ^,id,fc^{(s l+ i, Jd ,r I+1 , fc j}(^) + 9i > id > fe/{(s i+ i, Jd ,r- i+ljfcd )}(^), 

where <5{ a } denotes here the Dirac mass in a 6 I 2 and the above probabilities 
are given in Proposition El 
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We set now (S^, r£ )te[o,T] as the continuous-time process defined through the linear 
interpolation (in time) of the chain: for t G [ih, (i + l)h], 

^ = # + ^(*i-#) and f^ + ^^-r*). 

We are now ready to state the convergence in law of the Markov chain (rf, )i=o,...,jv 
to the diffusion process of our interest. 

Theorem 7 The Markov process (S^,f^) t e[o,T] converges in law on the space of the 
continuous functions C([0,T]; [0, +oo) x [0, +oo)) endowed with its Borel a-algebra 
to the diffusion process (S t , r t )te[o,T] solution to the equations ([[]) and ([2]). 



Proof. The proof is standard, see e.g. Nelson and Ramaswamy [10] or also classical 
books such as Billingsley [3], Ethier and Kurtz [7] or Stroock and Varadhan [TTj . 

To simplify the notations, let us set 

Mf dtk (£) = E((5* t - S } lY | (S* r*) = [S,,. ,;,)). i = 1,2,4, 
Ml hk {i) = E((r* +1 - r}Y | (Sf,r?) = (S id ,n, k )), t = 1,2,4, 
MfSjk = HiShi - - r?) | (S* r?) = (S,, r ,;.,)) . 

It is clear that Aifj k (£) is the local moment of order i at time ih related to S, as 

well as -M[ fc (f) is similar but related to the component r, and M-i'j >k is the local 
cross-moment of the two components at the generic time step i. 
So, the proof of the theorem relies in checking that for r* > and S 1 * > fixed, 
setting A* = k) : r i)k < r*, Sij < S*}, then the following properties i), ii) and 

in) hold: 



i) (convergence of the local drift) 

1 (ij',fc)eA» 



lim sup ^ |-M^- fc (l) - (/is)i,i,fe/i| 



lim sup - \M r i jk (l) - (/i r )i,fc/i| = 0; 

{i,j,k)£A, h 



ii) (convergence of the local diffusion coefficient) 

1 



Ii',',', -up - \Mt J)k (2) - a^hl = 0, 
1 

(iJ,fc)6A, 



lim sup - \Ml jk (2) - a r 2 r iifc /i| = 



lim sup r \Mi'h - p<T r <TsSi,j^/rXk~h\ = 0; 
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Hi) (fast convergence to of the fourth order local moments) 

lim sup -Mf Jtk (A) = 0, 
lim sup -Ml J)fc (4) = 0. 

So, we consider h such that h < min(/ii, /12, ^3), h±, hi and /13 being given in Lemma 
121 Lemma H] and Proposition [6] respectively, so that we can use both Remark [3] and 
Remark 

Proof of i). By using fl30|) and (137]) . we immediately get that -M^- jfe (l) = {[is)i,j,kh 
and A^[j fc (l) = (fi r )i,kh, and i) holds. 

Proof of ii). As for the cross-moment, this is really immediate when r^ k > 9*\/~h: the 
transition probabilities solve system ( 12T]) . whose last equation is actually ■M i 'j k = 
pcr r crsSij^/f\kh. Consider now the case rj fc < 9*\fh: here the transition probabilities 
are of the product-type, so that -Mf^, = (^ r )i,k{f^s)i,j,kh 2 ■ Therefore, 

\^i,j,k ~ P a rO-sSi,j^/r~^h\ = {^r)i,k{ps)i,j,kh 2 + o- r o- s Sijy/r~i~h 

< k(9 + r*)r*S*h 2 + arasS^^/^h 5 ^ 

< Ch 5 /\ 

As for the second order moment, again by using (130|) and (137]) . it follows that 

-^i,fc( 2 ) = (Si+i,ju + Si+i,id ~ 2Si,j)(f J 's)i,j,kh + (S i+1:ju - Sij)(Sij - S i+lijd ) 
M[ k (2) = (r i+1:ku + r i+hkd - 2r ijk )(n r ) ijk h + (r i+1>ku - r i>k ){r ijk - r i+1>kd ). 

So, by using (136|) . we get 

Ml hk {2) = Sl^r^cosHasVh) -l)h + (e^ - l) 2 ). 

Therefore, for (i, j, k) G A* we get 

\M? Jtk (2) - a 2 s Slh\ < 2S 2 M™M°sSh) -l)h + S 2 \e~° sV ~ h (e^ h - l) 2 - a 2 h 

and the statement holds. Concerning the second moment on r, we first study the 
case r,- k as in Lemma [21 Then by using ( 127]1 we have 

\Ml k (2) - a 2 .r i)k h\ < 2C*K0h 3 / i+1 + C 2 h 3/2 < Ch 3 / 2 . 

Consider now the case 0*Vh < r i>k < r*. Then, for h < (9*/r*) 2 then r ijk < 8*/y/h, 
9* as in Lemma [21 So, we use (12 8p and we obtain 

\MU2) - a 2 r hk h\ < k(9 + h 2 + ^h 2 
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and the statements again holds. 

Proof of in). First, straightforward computations give 

= ((^s)i,j,kh + Sij - S i+1:jd )(S i+ ij u + S i+ltjd - 2Sij)x 
x ((Si+i,j u — Sij) 2 + (S i+ x > j d — Sij) 2 ) + (S i+ ij d — Sij) A 
•Mi,fc( 4 ) = {{^ r )i,kh + r i)k - r i+ltkd )(r i+ltku + r i+hkd - 2r itk )x 

x ((r i+ i, fc „ - r i)k ) 2 + (r i+ i ;fed - r itk ) 2 ) + (r i+hkd - r iifc ) 4 . 

As for the first quantity, we immediately get 

Mfj >k (4) <4:St(nh + 1 - e~ asVK )(cosh(a s Vh) - l)2(e asVK - 1) 2 + 
+ St(e- asVK - l) 4 < Ch 2 

Concerning the 4th moment for r, first notice that if r^k < 9*vh then ( 1271) gives 
•M[ fc (4) < C*/;. 3 . If instead < < r* then from fl28|) one gets 

/ <T 2 \ / / O 2 \ 2 / / cr 2 \ 4 

A4[ fc (4) <[k(6 + r*)h + a r ^r*h + -jhja 2 r h[a r ^r if h + -jhj + ^o- r Vr*/i+^/iJ 

<Ch 2 . 

So, the proof is completed. 
□ 



6 Numerical results 

In this section we compare the performance of our lattice algorithm (called ACZ) 
with the procedures of Wei (WEI) and of Hilliard, Schwartz and Tucker (HST) for 
the computation of European and American options with the CIR stochastic interest 
rate. 

In the European and American option contracts we are dealing with, we consider a 
set of parameters already used in Hilliard, Schwartz and Tucker [8]: So = 100, 
a s = 0.25, r = 0.06, 9 = 0.1, k = 0.5, p = -0.25. In order to study the 
numerical robustness of the algorithms, we choose different values for o r : we set 
o r = 0.08,0.5,1,3. Let us remark that for a r = 0.5,1,3, the Novikov condition 
2k9 > a 2 is not satisfied at all. The parameters of the option contracts are the 
following: the strike is K = 100 and the maturity T is varying, since we set T = 1,2 
years. Finally, the number of time steps N varies: N = 50, 100, 150, 200, 300. 
Tables [T] and [2] report European put option prices for T = 1 and T = 2 respectively. 
The benchmark value is obtained using a Monte Carlo with very large number 
of Monte Carlo simulation (10 million iterations) using the accurate Alfonsi [1] 
discretization scheme for the CIR process with M = 300 discretization time steps 
(this method provides a Monte Carlo weak second-order scheme for the CIR process, 
without any restriction on its parameters). We also provide results for American 
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put option prices, as reported in Table [3] and Table H] (no benchmarks are available 
in this case). 

The numerical results show that our method provides very reliable and stable out- 
comes. Even if no odd results appear as as varies, on the contrary both WEI and 
HST fail when a r increases. As already observed by Tian [TJ], P3], this follows 
from the procedure set up to approximate the CIR process and can be explained 
by looking at the behavior of the drift (see 01])) associated to the transformed 
process R, which is the one used to define the transition probabilities (see Remark 
[[]). Indeed, one can write 

, \ 1 fal-AnO 

MR) = -tttH-^ — + Ra r 



2af V R 

As a r = 0.5, 1, 3, one gets of — An6 > 0, so that when R \, one has Hr(R) | — oo 
and pi,k — > 0. This gives that, in some sense, the process tends to be absorbed in 0, 
so that from the numerical point of view it becomes crucial to drastically decrease 
h. But this brings to procedures which are computationally heavy from a memory 
requirement point of view. 







N 


WEI 


HST 


ACZ 


MC Benchmark 




=0.08 


50 


6.599590 


6.579999 


6.547556 








100 


6.596199 


6.586245 


6.569844 


(6.580864) 






150 


6.594964 


6.588312 


6.577486 


6.586622 






200 


6.594367 


6.589255 


6.581016 


(6.592380) 






300 


6.593705 


6.590247 


6.584744 




0> 


=0.50 


50 


6.522023 


6.495957 


6.515875 








100 


6.606752 


nan 


6.543780 


(6.544551) 






150 


6.708564 


6.430882 


6.551733 


6.550315 






200 


6.588212 


6.308088 


6.556571 


(6.556078) 






300 


6.492437 


6.280427 


6.560239 






=1.00 


50 


7.397553 


3.640714 


7.054686 








100 


4.786837 


nan 


7.109525 


(7.153304) 






150 


8.046617 


4.094663 


7.123767 


7.159471 






200 


4.846171 


4.082429 


7.127271 


(7.165637) 






300 


7.340537 


4.135390 


7.126012 






=3.00 


50 


9.582836 


0.042281 


8.491568 








100 


3.026463 


0.074981 


8.636541 


(8.756826) 






150 


15.966217 


0.085990 


8.671644 


8.763625 






200 


9.103675 


0.081284 


8.648793 


(8.770423) 






300 


1.604996 


0.084082 


8.681638 





Table 1: European put options with T = 1, Sq = 100, K = 100, as = 0.25, tq = 0.06, 
9 = 0.1, k = 0.5, p = —0.25, a r varying. 
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N 


WEI 


HST 


ACZ 


MC Benchmark 


oy 


=0.08 


50 


7.113252 


7.075978 


7.044844 








100 


7.109863 


7.090808 


7.075287 


(7.090164) 






150 


7.108682 


7.095834 


7.085454 


7.096171 






200 


7.108057 


7.098328 


7.090541 


(7.102178) 






300 


7.107381 


7.100858 


7.095698 




(T r 


=0.50 


50 


7.590023 


nan 


7.532236 








100 


7.384775 


7.151360 


7.573197 


(7.575268) 






150 


8.099945 


7.157489 


7.590692 


7.581702 






200 


7.866820 


nan 


7.598707 


(7.588135) 






300 


8.627843 


7.029057 


7.600614 






=1.00 


50 


10.511291 


nan 


9.029287 








100 


10.377910 


1.513991 


9.113234 


(9.312621) 






150 


11.204261 


1.603589 


9.204665 


9.319903 






200 


3.979039 


nan 


9.192707 


(9.327185) 






300 


12.541620 


1.839346 


9.206984 




oy 


=3.00 


50 


0.003130 


0.000216 


11.801614 








100 


13.515823 


0.000344 


11.656249 


(12.046144) 






150 


12.799415 


0.000594 


11.820232 


12.054222 






200 


1.864718 


0.000832 


11.873238 


(12.062299) 






300 


24.377573 


0.001110 


11.919532 





Table 2: European put options with T = 2, Sq = 100, K = 100, as = 0.25, ro = 0.06, 
9 = 0.1, k = 0.5, p = —0.25, oy varying. 





N 


WEI 


HST 


ACZ 


a r =0.08 


50 


7.459623 


7.422693 


7.438208 




100 


7.456394 


7.437777 


7.445392 




150 


7.455096 


7.442827 


7.447790 




200 


7.454461 


7.445104 


7.448865 




300 


7.453750 


7.447485 


7.449971 


a r =0.50 


50 


7.647938 


7.596269 


7.650304 




100 


7.698584 


nan 


7.662723 




150 


7.751229 


7.586111 


7.665955 




200 


7.684290 


7.472087 


7.668003 




300 


7.620723 


7.627011 


7.669464 


a r =1.00 


50 


8.187062 


6.511775 


8.084019 




100 


6.750600 


nan 


8.106474 




150 


8.743109 


6.506826 


8.113036 




200 


6.822003 


6.501217 


8.115025 




300 


8.164598 


6.534435 


8.116760 


o r =3.00 


50 


9.699021 


2.902701 


8.953802 




100 


5.821256 


3.165914 


9.009359 




150 


16.069688 


3.405395 


9.028517 




200 


9.212406 


3.307452 


9.025299 




300 


4.740153 


3.049339 


9.037878 



Table 3: American put options with T = 1, Sq = 100, K = 100, as = 0.25, ro = 0.06, 
9 = 0.1, k = 0.5, p = —0.25, a r varying. 
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N 


WEI 


HST 


ACZ 




=0.08 


50 


9.161152 


9.087704 


9.149096 






100 


9.162292 


9.125087 


9.155979 






150 


9.162352 


9.137484 


9.158033 






200 


9.162318 


9.143667 


9.159035 






300 


9.162238 


9.149867 


9.160028 


(J r 


=0.50 


50 


9.814438 


nan 


9.821675 






100 


9.717844 


9.555458 


9.833291 






150 


10.113582 


9.570541 


9.838509 






200 


9.986565 


nan 


9.841491 






300 


10.358344 


9.576959 


9.842231 




=1.00 


50 


11.421349 


6.586697 


10.790323 






100 


11.365461 


6.934987 


10.834057 






150 


11.826422 


6.392858 


10.868361 






200 


7.790298 


nan 


10.863552 






300 


13.160854 


6.897618 


10.871842 




=3.00 


50 


2.837551 


3.291170 


12.488520 






100 


13.585301 


2.903662 


12.458445 






150 


12.865822 


2.672643 


12.516032 






200 


6.335779 


3.166955 


12.539396 






300 


24.420598 


3.405717 


12.564801 



Table 4: American put options with T = 2, S = 100, K = 10, a s = 0.25, r = 0.06, 
6 = 0.1, k = 0.5, p = —0.25, a r varying. 



21 



References 



[I] Alfonsi, A. (2010): High order discretization schemes for the CIR process: appli- 
cation to affine term structure and Heston models. Mathematics of Computation, 
79, 209-237. 

[2] Amin, K., Khanna, A. (1994): Convergence of American option values from 
discrete-to continuous-time financial models. Mathematical Finance, 4, 289-304. 

[3] Billingsley, P. (1968): Convergence of probability measures. Wiley, New York. 

[4] Brennan, M. J., Schwartz, E. S. (1976): The pricing of equity-linked life insur- 
ance policies with an asset value guarantee. Journal of Financial Economics, 3, 
195-213. 

[5] Briani, M., Caramellino, L., Zanette, A. (2013): Pricing American options in 
the Heston model: an hybrid tree-finite difference approach. Working paper. 

[6] Cox, J.C., Ingersoll, J., Ross, S. (1985): A theory of the term structure of interest 
rates. Econometrica, 53, 385-407. 

[7] Ethier, S.N., Kurtz, T. (1986): Markov processes: characterization and conver- 
gence. John Wiley & Sons, New York. 

[8] Hilliard, J.E., Schwartz, A.L., Tucker, A.L. (1996): Bivariate binomial pricing 
with generalized interest rate processes. The Journal of Financial Research, XIX- 
4, 585-602. 

[9] Kushner, H., Dupuis, P.G. (1992): Numerical methods for stochastic control 
problems in continous time. Springer Verlag. 

[10] Nelson, D. B., Ramaswamy, K. (1990): Simple binomial processes as diffusion 
approximations in financial models. The Review of Financial Studies, 3, 393-430. 

[II] Stroock, D. W., Varadhan, S.R.S. (1979): Multidimensional Diffusion Pro- 
cesses. Springer, Berlin. 

[12] Tian, Y. (1992): A simplified binomial approach to the pricing of interest-rate 
contingent claims. Journal of Financial Engineering, 1, 4-37. 

[13] Tian, Y. (1994): A reexamination of lattice procedures for interest rate- 
contingent claims. Advances in Futures and Options Research, 7, 87-110. 

[14] Wei, J., Z. (1996): Valuing American equity options with a stochastic interest 
rate: a note. The Journal of Financial Engineering, 2, 195-206. 



22 



